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Hydrodynamic flow in both classical and quan- 
tum fluids can be either laminar or turbulent. 
To describe the latter, vortices in turbulent flow 
are modelled with stable vortex fllaments. While 
this is an idealization in classical fluids, vortices 
are real topologically stable quantized objects in 
superfluids. Thus superfluid turbulence [1] is 
thought to hold the key to new understanding 
on turbulence in general. The fermion super- 
fluid ■^He ofl"ers further possibilities owing to a 
large variation in its hydrodynamic characteris- 
tics over the experimentally accessible tempera- 
tures. While studying the hydrodynamics of the 
B phase of superfluid ■^He, we discovered a sharp 
transition at 0.60 Tc between two regimes, with 
regular behaviour at high-temperatures and tur- 
bulence at low-temperatures. Unlike in classi- 
cal fluids, this transition is insensitive to velocity 
and occurs at a temperature where the dissipa- 
tive vortex damping drops below a critical limit. 
This discovery resolves the conflict between exist- 
ing high- and low-temperature measurements in 
■^He-B: At high temperatures in rotating flow a 
vortex loop injected into superflow has been ob- 
served to expand monotonically to a single rec- 
tilinear vortex line [2 4], while at very low tem- 
peratures a tangled network of quantized vortex 
lines can be generated in a quiescent bath with 
a vibrating wire [5]. The solution of this conflict 
reveals a new intrinsic criterion for the existence 
of superfluid turbulence. 

In conventional liquids the vorticity u; = V x v obeys 
the Navier-Stokes equation 



— V X [v X u;] 



conditions through the Reynolds number, Re = RU/v, 
formed by the characteristic size R of the system, the 
flow velocity U, and the kinematic viscosity f. At large 
Re 3> 1, the effect of the inertial term is dominating, and 
laminar flow becomes increasingly unstable. If vorticity is 
released into a meta-stable laminar state, a sudden tran- 
sition to a chaotic flow of eddies occurs. The evolution of 
the turbulent flow is described by the Kolmogorov energy 
cascade: the kinetic energy of the flow is transferred to 
ever smaller length scales, with large vortex loops decay- 
ing into smaller loops, until a scale is reached where the 
energy is dissipated by viscosity. 

In superfluids turbulence acquires new features: A 
superfluid can be described as consisting of two inter- 
penetrating fractions, the frictionless (superfluid) and 
viscous (normal) components. If both fractions are mov- 
ing, as is the case eg. behind a pulled grid in measure- 
ments on superfluid turbulence in ^He [1], the turbulent 
state bears more resemblance to that of viscous liquids. 
A new class of superfluid turbulence becomes possible 
when the normal component is so viscous that it is es- 
sentially immobile and only the superfluid component is 
moving with respect to the boundaries. This is the usual 
case in the flow of "^He, which is considered here. 

The vorticity of the superfluid component is quantized 
in terms of the circulation quantum k = 2Trh/M where 
M is the mass of a superfluid particle. The flow of the 
superfluid component can be characterized with the 'su- 
perfluid Reynolds number' Rcs = RU/k. The Feynman 
criterion Rcg ~ 1 gives the velocity at which it becomes 
energetically favorable to form a quantized vortex. If a 
large nucleation barrier exists, vortices are not created 
and the superfluid remains vortex-free, in our measure- 
ments here up to Rcs ^ 200. Quantized vortices can then 
be injected, either by some extrinsic means or by increas- 
ing the velocity so that the nucleation barrier is overcome 
and spontaneous vortex formation occurs. The limit of 
large Rcg 3> 1 is equivalent to a vanishing Planck's con- 
stant, K oc /i — > 0, so that the vorticity becomes a con- 
tinuous variable, like in the classical case. However, the 
fluid still remains unconventional because of its two-fluid 
nature: In addition to the superflow, there is the nor- 
mal component which is at rest in the container frame, 
Vji = 0. Interactions between the superfluid vorticity and 
the normal component give rise to a mutual friction force 
on a unit volume of the superfluid. 
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The interplay of the two terms on the rhs, the inertial 
first term and the viscous second term, governs the tran- 
sition to turbulence [6]. It is determined by the external 



which is described by the dimensionless temperature- 
dependent parameters a and a'. These correspond to 
the usual mutual friction parameters, if we are dealing 
with one vortex or locally polarized vorticity, or they 
may be renormalized in the general case. Here Vg and 
= V x Vg are the superfluid velocity and vorticity, 
and a>s is the unit vector in the direction of LJg- Insert- 
ing the mutual-friction force in the Eulcr equation one 
obtains the hydrodynamic equation for Rcg 3> 1 [7]: 
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FIG. 1. Summary of events: A vortex half ring is injected 
into vortex-free supcrfiow, generated with rotation around the 
vertical axis (f2 = 0.21rad/s or Res « 30). The figures rep- 
resent snapshots from simulation calculations in the rotat- 
ing frame, which are consistent with the experimental ob- 
servations. The only temperature dependence is contained 
in the mutual-friction parameters a and a' which wc take 
from Ref. [10]. (Top) At high temperatures the loop grows 
monotonically into a single rectilinear vortex line, pulled by 
the Magnus force from the superflow. The motion is highly 
damped since the Magrms force is balanced by the mutual 
friction force. (Bottom) At low temperatures the loop moves 
with the circulating superflow: (44 s) The vortex is unstable 
with respect to loop formation and develops through recon- 
nections into a tangle at the injection site. (68 s) The Magnus 
force from the vortex-free superflow extends the tangle along 
the vertical axis. The measured flight time tf for the tangle 
to extend across a vertical distance d is controlled by mu- 
tual-friction damping: tf ~ d/{Q.Ra). (108 — 160s) The 
final stop is the roconnoction of the tangle to rectilinear lines. 
The growth in the relative polarization is experimentally ob- 
served as a relaxation in the NMR absorption and the total 
number of vortex lines in the final stable state is determined 
from the NMR Une shape (Fig. 2). 



dt 



= (1 - a') V X [vs X a>s] + aV x [Cj^ x (a;^ x vj] . 

(1) 



The energy dissipation is determined by the mutual fric- 
tion damping a in the second term, while the reactive 
mutual-friction a' renormalizes the inertial term of con- 
ventional hydrodynamics. Like in a classical fluid, the 
inertial term drives the instability by increasing the num- 
ber of vortex loops, whereas the dissipative term acts to 
stabilize laminar flow by decreasing the number of loops. 
The fundamental difference from conventional hydrody- 
namics is that these two competing terms now have a sim- 
ilar dependence on the velocity and its gradients. Their 



relative importance is determined by the intrinsic dimen- 
sionlcss parameter of the superfluid, q = 0^/(1 — oi'), in 
contrast to classical liquids where it is governed via the 
Reynolds number by the extrinsic quantities R and U. 
We thus expect laminar flow to become unstable when q 
drops below a critical value qc ^ \- In the case of a single 
vortex line the parameter q also marks the crossover from 
propagating Kelvin waves [q < 1) to the over-damped 
regime {q^ 1) [8]. In Fermi supcrfluids and superconduc- 
tors with finite energy gap one has q « (wqt)"^, where ujq 
is the spacing between the bound states of quasiparticles 
in the vortex core and is their lifetime broadening, 
owing to scattering from the normal component [9]. In 
^He-B, q approaches infinity at and drops monotoni- 
cally to zero, when cooled down to T — > [9,10]. Such an 
evolution of q with temperature would in classical fluids 
correspond to scanning the Reynolds number from to 
oo: We can thus study the emergence of turbulence as a 
function of temperature when q falls below unity. 

Measurements: — In rotating "^He-B high-velocity 
vortex-free superflow can be achieved, into which vor- 
tex seed loops can be injected (Fig. 1). We record as a 
function of temperature T and rotation the number 
of rectilinear vortex lines after injection (A^f), when all 
transients have died out. Non-invasive NMR measure- 
ment is used to determine the number of vortex lines in 
the final state. The NMR absorption spectrum of ^He-B 
maps the order-parameter texture, which is strongly af- 
fected by the orienting effect from the Vg field. It allows 
us to determine TVf with a resolution better than 10 vor- 
tex lines in the temperature range of Fig. 2. In Fig. 2d 
we list the result after each injection event in the (fi, T) 
plane, classified as a regular high-temperature (blue) or 
turbulent low-temperature (red) process. Here the in- 
jection method is the shear-fiow instability of the phase 
boundary between the ^He-A and ^Hc-B phases, which 
produces a small random number (AA^ ~ 10) of B-phase 
vortex loops whose one end sticks out of the AB interface 
while the other end is on the cylindrical sample boundary 
[4]. At high temperatures an injected loop expands and 
becomes a single rectilinear vortex line. At low tempera- 
tures the loop evolves into tangled vorticity which finally 
fills the whole sample volume with rectilinear lines. The 
central plot Fig. 2d leads to a surprising conclusion: Only 
a narrow cross-over region at 0.60 Tc separates the two 
processes, with little dependence on the maximum su- 
perflow velocity U = ClR of the initial vortex-free state. 
At the transition temperature 0.60 the experimental 
value of g = a/(l - a') « 1.3 [10]. 

Vortices can also be injected with other methods. In 
the absence of the AB interface vortices are created (i) at 
a rough spot on the cylindrical sample boundary [2] , (ii) 
under neutron irradiation [3] , or (iii) by a sudden leak of 
vortex loops through the small oriflce which connects the 
sample to the rest of the liquid ^He volume [11]. These 
three processes have different critical velocities and tem- 
perature dependences, but in all cases the measurements 
are consistent with a crossover at 0.60 Tc from a low-yield 
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FIG. 2. Principle of measurement and central result: a, The initial state is vortex-free superflow in rotation, with the normal 
component stationary and the superfluid component moving in the rotating frame. This state is identified by a large sharp 
peak in the NMR absorption spectrum, which is caused by the vortex-free superflow and is shifted far from the Larmor value 
(line shape denoted as M = in b and c). When a few (AA'') vortex loops are injected, the number of vortex lines in the final 
steady state (A^f) is found to fall in one of two categories: b, A^f = A7V, indicated by a small shift of the NMR absorption 
from the superflow peak close to the Larmor value. The shifted absorption is proportional to Nf, as shown in the insert by 
the build up of absorption at small frequency shift. This process corresponds to regular loop expansion, c, AA'^ Nf < A^cq, 
indicated by a total removal of the superflow peak (red). Here A'eq ~ 7r7?^2f2/K ~ 10^ is the equilibrium number of vortex 
hnes in rotation. This process corresponds to turbulent loop expansion. The relaxation of the line shape through intermediate 
stages (c) represents a direct signal from the turbulent epoch: In the final state with a large number of rectilinear vortex lines 
{N{ ~ A'eq), the distribution of absorption is relatively flat. In contrast, in the turbulent state with unpolarized Vs field, much 
of the absorption is piled in a narrow peak at the Larmor edge. The intermediate curves between these two limiting line shapes 
represent the transient relaxation when the turbulent vorticity with uia ^ Weq decays to a polarized cluster of rectilinear lines 
with oJcq ~ 257. d, Phase diagramme of measured events: Each data point monitors the outcome from an injection event which 
is initiated using the shear-flow instability of the phase boundary between '^He-A and '^He-B. Vortex loop expansion separates 
into regular events (blue) and turbulent events (red) as a function of rotation Q, and temperature T, with an abrupt transition 
at 0.60 Tc independently of Q. Only within a very narrow interval AT ~ 0.05 Tc around 0.60 Tc can intermediate values of TVf 
be measured. By scanning the temperature (curved lines) or the current in the magnet for stabilizing the AB phase boundary 
(vertical lines), the critical velocity Vc ~ of the interface instability can be varied along well-defined continuous trajectories. 



to a high-yield process with decreasing temperature. De- 
pending on sample geometry, below 0.50 Tc the vortex 
leak through the orifice may occur at 0.5rad/s or less. 
Even at these low velocities {Res < 70) vortex loop ex- 
pansion is found to be turbulent. 

The long cylindrical sample is monitored with two sig- 
nal coils at different locations, connected to indepen- 
dently operated NMR spectrometers. In addition to the 
number of vortex lines in the final stable state, we can 
obtain information on the transient state: (i) on the axial 
expansion of the vorticity and (ii) below 0.60 Tc on the 
decay of the tangle into a polarized array of rectilinear 
lines. The former process, the expansion of the vorticity 
along the tt axis between the two signal coils, is con- 



trolled by the mutual-friction damping a: The extracted 
flight time for the motion in the axial direction repro- 
duces a{T) from Ref. [10] and is continuous across the 
transition at 0.60 Tc. Thus the axial expansion is indif- 
ferent to whether the vorticity moves as nonintcracting 
loops or as a network of loops. 

Below 0.60 Tc, the arrival of the tangled vorticity 
within a pick-up coil is signalled by a massive abrupt 
shift of the absorption close to the Larmor edge, where 
a new sharp absorption peak is formed (Fig. 2c). This 
line shape is similar to that measured in the non-rotating 
state (0 = 0). Here the orientational effects from the su- 
perflow fields trapped around the vortex cores are aver- 
aged out, owing to the randomly distributed vortex tan- 
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glc with density Us ^ t^^cq and the removal of essentially 
all large-scale vortex-free superflow. The sharp peak then 
decays with a relaxation time ~ 30 s towards the broad 
final stable state spectrum of a polarized vortex array 
with rectilinear lines and Ug = Weq = 2f2. This restruc- 
turing of the NMR lino shape is thus a direct signal from 
the transient epoch with a vortex tangle. 

Numerical simulation — As seen in Fig. 1, the measured 
features are consistent with our simulation calculations. 
We use the vortex filament model [12] in the rotating 
frame [13]. A vortex is represented in parametric form 
by s = s(^,t), where s refers to a point on the filament 
and ^ is the arc length along it. The spatial and time 
evolutions are integrated rigorously using the Biot-Savart 
law. The vortex velocity s is calculated from the dynamic 
equation s = Vgi -I- a's' x [s' x Vgi] —a[s'x Vgi] (recall that 
Vn = 0), where the local superfluid velocity Vgi includes 
all contributions to the superflow at s(^, t). As boundary 
conditions we use smooth solid- walls with image vortices. 

The initial state is a vortex half ring (diameter 3.5 mm) 
in the centre of a cubic sample container (10 mm wide). 
At high temperatures T/Tc = 0.8 the high damping a 
causes the seed loop to expand, mostly by motion in the 
radial and axial directions, into a single rectilinear vor- 
tex line in the centre of the sample. At low temperatures 
T/Tc = 0.4 the seed loop travels azimuthally and under- 
goes the Kelvin-wave instability, when the superflow is 
oriented along the vortex [8] . Here g < 1 so that Kelvin 
wave oscillations are only lightly damped and loop for- 
mation starts via reconnect ions. Thus a vortex network 
develops already at the injection site. The vortcx-frec 
superflow above the tangle stretches the loops with right 
orientation, while others shrink, and thus the tangle ex- 
tends and travels in the axial direction. Through recon- 
nections the properly oriented loops are gradually joined 
to lines which form a bundle almost parallel to the ro- 
tation axis Jl. The bundle precesses initially around the 
Ct axis, but presumably more reconnections of remaining 
loops will continue adding lines to the bundle and thereby 
reducing the superflow around the bundle, until the equi- 
librium number of rectilinear vortex lines is reached. The 
complete calculation is so far too time consuming. 

Discussion: — Both the results from the numerical sim- 
ulations in Fig. 1 and from the measurements in Fig. 2 
support the conclusion from Eq. (1) that the transi- 
tion to turbulence is determined by the parameter q = 
a /{I — a'): There exists a critical value gc ~ 1 which 
separates the mutual-friction-dominated regime a,t q > 
from the inertia-dominated turbulent regime at q < q^. 
With q> qc, the vortex- loop expansion is a regular stable 
process even in high- velocity superflow with Rcg ^ 1, as 
has been verified in many previous measurements. With 
q < qc, & multiplication process is switched on, the vor- 
tex lines create loops, become entangled, and reconnect 
forming new loops. Could qc be a fundamental universal 
number which, when ^ 1, places the upper bound 
to superfluid turbulence? 

Why has such a transition, driven by an intrinsic pa- 



rameter, not been observed previously in superfluid ^He? 
A number of reasons can be listed: (i) In ■'He the viscosity 
of the normal fraction is four orders of magnitude higher. 
Thus there is no ambiguity about the state of the normal 
component, unlike in the case of '*He. (ii) The vortex- 
core radius in ^He-B exceeds by two orders of magnitude 
that of ''He, where it is only of atomic size. Therefore 
vortex pinning and remanence at solid walls can often 
be neglected in "^Hc-B and metastablc states with rapid 
superflow become possible, (iii) The most important dif- 
ference lies in the value of mutual friction. In ^He the 
dissipativc part a is small and influences vortex motion 
relatively little. Only within a narrow temperature in- 
terval AT/Tx ~ 2 • 10-3 from Tx one finds g > 1 and 
expects different behavior. No conclusive measurements 
exist from so close to Tx- 

In the helium superfluids, the transition at ~ 1 lies 
very close to Tx in "^He, in the middle of the experimen- 
tally accessible temperature range in ''He-B, and below 
that at extremely low temperatures in "^He-A. This pre- 
diction is consistent with present experimental informa- 
tion. Here we found that turbulence is unstable even at 
high Rcs, if q > qc- Perhaps in other miflti-componcnt 
hydrodynamic systems, superfluid or viscous, the stabil- 
ity of turbulence might also be governed by an intrinsic 
parameter, in addition to the Reynolds number. 
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